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Abstract. We provide cubature formulas for the calculation of derivatives of 
expected values in the spririt of Terry Lyons and Nicolas Victoir. In financial 
mathematics derivatives of option prices with respect to initial values, so called 
Greeks, are of particular importance as hedging parameters. The proof of 
existence of Cubature formulas for Greeks is based on an argument, which leads 
to the calculation of Greeks in an asymptotic sense — even without Hormander's 
condition. Cubature formulas then allow to calculate these quantities very 
quickly. Simple examples are added to the theoretical exposition. 



1. Introduction 

Cubature formulas provide approximative values for integrals with respect to a 
given measure. The well-developed theory of cubature formulas in finite dimensions 
was recently applied to provide Cubature formulas on Wiener space by Terry Lyons 
and Nicolas Victoir in the beautiful seminal article 

We try to extend the framework of Terry Lyons and Nicolas Victoir to the cal- 
culation of Greeks. We briefly outline in the introduction the main results and 
possible applications of it. In the sequel we shall always work with C°°-boundcd 
vector fields Vi and C°°-bounded function /. Given a stochastic differential equa- 
tion in K N of the type 

d 

dY t y = V (V?)dt + ]T ° dBt, 

i=i 



on a stochastic basis (fi, (Ft)o<t<T, P) with d-dimensional Brownian motion, then 
cubature formulas provide a method to approximate E(f(Y^)) 1 namely - choosing 
a degree m of accuracy - there is number r > 1 and there are 7J 1 -trajectories 
ojj : [0, T] -> R d+1 and weights Ay > 0, j = 1, . . . , r, such that 

r 

E(f(Y?)) = £ Ai/Wfo)) + 0(t^), 

3=1 
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where Y^(ojj) denotes the point in M , which is obtained by solving the ordinary, 
non-autonomous differential equation 

i=0 

The constant in the order estimate depends in general on derivatives up to order 
m + 1, in the hypo-ellitptic case this dependence can be improved to the first 
derivative. Note also the important fact that the number r of trajectories can be 
bounded by a number depending on d and to, but not on ./V (see [8] for the precise 
estimate). Obviously this procedure can be iterated with small time steps by the 
Markov property, which yields then a high order numerical approximation scheme 
for the calculation of expected values E(f(Y^)). 

In the hypo-elliptic case we can calculate the derivatives of y i— > E(f(Y^)) with 
respect to the initial value by weight formulas, i.e. for any direction v £ R w there 
is a random variable 7r such that 

(1-1) ±U =0 E(f(Y* + ™)))=E(f(Y t V)ir) 

for bounded measurable /. Explicit formulas for ir can be given by Malliavin 
Calculus (see for instance jS] and and are in principle well-known since long 
time (see for instance [2]). 

If we fix an order of approximation m > 1 and a (homogenous) direction v, such 
that (for the notations see Section [5J 

v= M---[V il ,V i2 } > V i3 ],...,V ik ](y) 
ie.A\{0,(o)} 

dcg(/) = fc 

holds with some 1 < k < to, then we can construct a universal weight tt™ 1 (calcu- 
lated in a universal way from vi and the Brownian motions, see Remark [SJ, such 
that 

l\ e=0 E(f(Y t y +tv )))^ - E(f(Y?)-Kf tl )+0{t^) 

holds true, where the constant in the order estimate only depends on the first 
derivative of /. Theorem [S] asserts that we are able to find weights fij ^ and 
ff^trajectories Uj : [0, T] -> R d+1 such that 

(1-2) ^l«=oS(/(i? +eu ))** =^/i i /(y t B (^)) + 0(t^) 

holds true. Also here the constant in the order estimate depends on derivatives of 
/ up to order to + 1. The number r can be estimated similarly as in [S]. 

A combination of the original Cubature formula and formula H.2JI then yields 
by iteration a high order numerical approximation scheme, i.e. 

(1-3) j e \^ E{f(Yt ev )) = jl= E{E{f{Yr ev )\Fs )) 

(1-4) =±\^E{E{f{Y:_J)\ z=Y y^). 

Calculation of the inner expectation by the original cubature formula and of the 
outer expectation by formula H1.2JI lead to the desired iterative procedure: in the 
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uniformly hypo-elliptic case (see [B] for details) we consequently obtain the following 
recipe. 

• choose a small sq > and relatively large t — sq. Subdivide s = to < t\ < 
• • • < ffc = t with Si :— ti — for i = 1, . . . , k. 

• calculate by (|1.2|) with the function z i— > E(f(Y t f SQ )) for a certain degree of 
accuracy m. The order estimate yields (by careful choice of homogenous) 
directions 0{s m± ^~ jL ) for some 1 < k < m. The constant in the order 
estimate depends on the (m + l)st derivative of z i— ► E(f(Y t z _ s )), which 
by methods of Malliavin Calculus - can be reduced to dependence on the 
first derivative of /, i.e. there is an absolute bound given by 



since t — sq is bounded away from 0. 
• calculate z i— > E(f(Y t z _ s )) by Cubature formulas with certain degree of 
accuracy m' (see 8J. The order estimate, which again - by methods from 
Malliavin Calculus (see jH]) - can be reduced to dependence on the first 
derivative of /, i.e. there is an absolute bound given by 



• finally we obtain a high order approximation scheme for -^\ (i= oE(f(Y^ +ev )) 
through solving ordinary differential equations with estimate of the absolute 
error through 



The procedure only involves ordinary differential equations and knowledge 
on the cubature trajectories. By choosing m big enough, one can bound 
the difference m — k from below. 

Remark 1. A direct application of formula is often very slow numerically 

in the purely hypo-elliptic case (see |o]j- In the elliptic case there are very efficient 
ways to simulate it (see 0]^. Hence in particular the purely hypo-elliptic case con- 
stitutes an interesting field of applications, since one does not need to worry about 
invertibility of the covariance matrix. 

Remark 2. For all cubature formulas there are algebraic methods to determine the 
trajectories ujj and the weights Xj, or fj,j, respectively. These methods are based on 
calculus in nilpotent Lie groups, which constitute the local model for hypo-elliptic 
diffusions in the sense of Gromov (see for instance the very well written monograph 
[TT] for a general outline, and [7] for applications in the theory of - stochastic 
- differential equations). 



Cubature formulas on Wiener space (see the seminal article [S]) rely on the 
analysis of hypo-elliptic diffusions on free, nilpotent groups on the one hand (see the 
recent book pQ), on the other hand on stochastic Taylor expansion (see for instance 
[Hj). Furthermore the proof of a Cubature formula on Wiener space appears as an 
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application of Tchakaloff's Theorem on finite dimensional cubature formulas (see 
\V2\). We state these three main ingredients and show - in order to explain the 
mathematical background of Cubature formulas - a sketch of the proof (following 

El)- 

We fix a probability space (fi, (J-t)o<t<T, P) together with a d-dimensional Brow- 
nian motion (-BJ)o<t<T,i=i,...,d' We need the following notations for convenience: 

(1) We abbreviate by A the set of all finite sequences I := {%%,..., ik) € 
{0, . . . , d} fc , fc>0and we define a degree function on A by 

deg(ii, ...,i k ):=k + #{l with i t = 0}, 

which simply means that appearing 0s are counted twice. Additionally, we 
define dcg(0) = 0. We define a semi-group structure on A via 

(it, ...,ifc) * 0'i, . . := {ii, . . -,ik,h, ■ ■ -,3i), 

* (ii, ...,i k ) = (k, ...,ik) * = (*i, ■ ■ .,ifc)< 

(2) We denote by A™ 1 the free, nilpotent algebra with d+1 generators eo, . . . , e<j, 
i.e. the set of all non-commutative polynomials in those variables, such that 
the following nilpotency relations hold: if deg(7) > m, then e^e^ ■ ■ ■ &i k = 
for all I € A. Hence is a finite dimensional, non-commutative, 
real algebra with unit element 1 and we are given a grading via the de- 
gree functions, i.e. a monomial e^e^ ■■■&i k is said to have degree n if 
deg(«i, . . . , ifc) = n. Denote by W n the linear span of all monomials of 
degree n, then we obtain 

a& = e™ w n , 

furthermore W p W q C W p+q (where we define W p = for p > m due to the 
given relations), Wo = K • 1, so A 7 ] l 1 is a graded algebra. We denote the 
canonical projections of x on the subspaces W n by x n for n > 0. 

(3) On the finite dimensional algebra A™ x we define the exponential series 

oo ^ 

exp(z) -XItT' 

where the series converges everywhere due to the nilpotency relations. We 
define the logarithm on elements x with Xq ^ 0. We identify xq with a real 
number, so the series 

log(x) = log(x ) + £ L¥±(£Z£0 )<. 

is well-defined, since x ~^° € ®™ =1 W„ (i.e. the series is finite). 

(4) We define the Lie algebra generated by eo, . . . , e<j with respect to the 
Lie bracket [x, y] := xy — yx. The Lie algebra inherits the grading from the 
algebra A™ 1 via U n := 0™ x H W n , for n > 0. Hence 

B™1 = ®™=l^ri 



is a graded Lie algebra. In fact the Lie algebra is free, nilpotent of step to, 
with d generators of degree 1 and one generator of degree 2. 
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(5) We denote the exponential image of q™ 1 by G^\ := exp(g™ 1 ) C A™ x and 
call it the free, nilpotent Lie group. G™ 1 is indeed a Lie group as closed 
subgroup of Lie group 1 + (B™ = iW n . The tangent space at x € G™ is 
spanned by the left (or right) translations xw for w 6 

(6) We need the canonical dilatations on Lie algebra and Lie group. We define 
for x S A™ x an algebra homomorphism A t via 

A t (x n ) = t n x n 

for x n € W n , n > 1 and t > 0. The homomorphism is well defined and 
restricts to a Lie algebra homomorphism A t on q™ 1 and and a Lie group 
homomorphism A t for t > 0. 

(7) We apply the notation odB^ := dt. We abbreviate the left invariant vector 
field associated to by D t , i.e. Di(x) — xei for x € G™ r We define a 
stochastic process (Xf )o<t<T on G^\ via 

(2.1) dX? = £A(X?)°dB* 



(2.2) =^Xf ei o^, 

i=0 

(2.3) = 

and see immediately that Xf — xX\ for < t < T almost surely. 
(8) Note here and in the sequel that vector fields on a vector space V : M. N — > 
M. N are used with double meaning: either as tangent directions on a smooth 
geometric object, or as first order differential operators on smooth functions 
/ : M w -> K, 

:=df(x)-^) 

for a; € R^. For smooth maps G : K w — > R M we apply the notion of the 
tangent map or Jacobian 

(2.4) dG(a:) •« := -^-| e=0 G(x + ew) 

ae 

for R . A vector field is called C 00 - bounded if all derivatives of order 

greater of order greater than are bounded. 

In order to see the relation between free, nilpotent Lie groups on the one hand 
and asymptotic analysis as t J, on the other hand, we formulate stochastic Tay- 
lor expansion (see pp and [3]), which simply results from iterating the defining 
equations for a stochastic differential equation. 

Theorem 1 (stochastic Taylor expansion). Given C 00 -bounded vector fields 

V ,...,V d :R N ^R N 
on M. N , then the diffusion process 

ii 



dY t y= J2 Vi(Y«)odBi, 



i=0 



yV - y 
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for y £ M. N admits the following series expansion: for n > and f £ C£°(R ) the 
expansion 

f{Y?) = £ V ix ■ ■ ■ V lk f(y) f odBH ° ' ' ' ° dB Z + 

T C A JO<t!<-<t k <t 



ieA 

dcg(/)< 



Rn(t,f)= J2 I V i0 ---V ik f(Y t y )odBi° o.-.o dBll 

,. \ — * J0<tn<---<t k <t 



deg(ii,...,i k )<n 

(i ,...,i k )&A 
deg(io, —,i k )>n 



holds true for < t < T. Additionally \J E(R n (t, f) 2 ) — 0(t 2 ) } where the 
constant in the order estimate depends on the (m + l)st derivative of f. 

Remark 3. The above procedure works for "any" differential equation (see the 
literature on rough paths, for instance the seminal work \J\), in particular, if we 
consider the non- autonomous, ordinary differential equation 

d 

dZ y t = ^V l {Z y t )duj\t) 

i=0 

for a ^-trajectory uj : [0,T] -> R d+1 , we have for n > and f £ C£°(M. N ) the 
expansion 

f(Zf) = V V h ■ ■ ■ V ik f{v) I dcj^ (h) ■ ■ ■ fa** (t k )+ 

IeA J 0<tx<-<t k <t 

dcg(/)<n 

+ Rn(tJ), 

Ru{t,f)= J2 I v i0 ---v tk f(zl)d^>(t Q )---du J ^(t k ), 

{i:,...,i k )C.A J0<t <-<t k <t 
dcg(ii,...,i fc )<n 

(i ,...,i k )eA 
dog(i ,...,ifc)>n 

holds true for < t < T. Notice here that in order to obtain an asymptotic 
expansion (for a given fixed trajectory u)) of certain order in time t, one has to 
change the degree function. 

Example 1. If we apply this series expansion to the process Xf in the vector space 
A™ 1; we obtain for any linear function A : A™ x — * R 



Apr?) - A ^ • ' ' e ^ ) / odB $l 

IeA Jo<t 1 <-<t h <t 

dcg(/)<m 

for x £ A^j and < t < T, since 



DiX(x) := — | e=0 A(x + exei) = X(xei) 
ae 

for a £ tf£ v hence 

(D h ■■■D ik X)(x) = \(xe tl ■■■e ik ). 
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Consequently we obtain in A™ x 

Xj= I odBl\o...odB\l 



leA 

dcg(/)<m 



0<ti<--<t fc <t 



for < t < T , almost surely, by duality. This formula provides a nice stochastic 
representation of the solution of eauation \2.1l where we apply that G™ x is embedded 
in the algebra A™ 1 . 

Stochastic differential equations provide solutions for the associated heat equa- 
tion via expectations, so we obtain for / g C^°(A^ l 1 ), 

(Do + lit -> E(f(Xm = §- t E(f(Xn), 

i—l 

for all x <E A j* 1 . The heat equation itself admits a similar construction as stochastic 
Taylor expansion, namely any solution u(t,x) with initial value u(0, x) = f(x), 
where / is C°°-bounded, can be written as 
n ■ d 

«(*• I )=E#»+oE + Rn(t, /), 



i 2 



/ • • ' / 1 (A) + \V\ D*) n+1 u(s , x)ds Q ...ds n . 
Jo Jo l l=1 

In the case of / = A for a linear function A : A™j — > M, this leads for the process 
in A™ 17 to the nice formulas 

£(A(Xf))= ^ X(xe il ---e ik )E([ odBfc o • • • o dfl£), 

j e _4 "'0<ti<-<t fc <t 
dcg(/)<m 

1 - 

u(t, x) = A(x exp(i(e + - ^E e i )))) 



2 

i=l 



and hence 

(2.5) exp(t(e + lj2 e "))= E e *i'" e i^(/ odB£o...odB&), 

i=l Je.A J0<ti<-<i fc <t 

dcg(/)<m 

for < t <T, which is one basic formula for the construction of Cubature formulas 
(see [I] and 0). 

In order to obtain a cubature formula we need a version of Tchakaloff 's theorem 
(see ^2] for details in the case where the support of /i is not compact): 

Theorem 2. Let n be a measure on M. N such that moments of all orders exist. Fix 
a number m, then there are points x\,...,x r G supp(/Lt) and weights Xj > for 
j = 1, . . . , r, such that 



I }{yMdy)=Y.tt x i) X i 
Js, N ~[ 

up to order m. Furthermore 
Po\ m (M. N ) is the vector space of polynomials on M. N with degree less or equal m. 



for all polynomials on M. N up to order m. Furthermore r < dimg Pol m (R Ar ), where 
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On basis of these prerequisites we can define cubature formulas on Wiener space: 

Definition 1. Fix m > 1 and < t < T, a cubature formula on Wiener space 
is given by a finite number of points x%,...,x r € G™ x and finitely many weights 
Ai, . . . , A r > 0, such that 

r 

E{X}) =Y,h*i> 

i=i 

or equivalently due to formula $2.5\) 

^ d r 

exp(t(e + - Y e i)) = Y X j x j- 
i=i j=i 

Furthermore r < diniR . 

Applying Tchakaloff's theorem to the law of the process X\ in Aj\, < t < T, 
which is supported in G™ 1; yields that we find points x\, . . . ,x r € supp((X t )*P) C 
G™ 1 and weights X 3 > 0, such that 

E{Xl)= Y, e n ---e lk E(f °dB%o ■ ■ ■ o dB^) =^X J x j , 

IeA Jo<ti<-<t h <t j=1 

dcg(/)<m 

whence the existence of Cubature formulas for any m > 1. 

Finally we apply Chow's theorem of sub-riemannian geometry (see pQ and 
which tells that every point x <E G™ x can be reached by a if -^horizontal curve, i.e. 
for every x £ G 1 J l 1 we find a ii^-curve ui : [Q,T] — > R d+1 , such that the solution of 
the non-autonomous ordinary differential equation 

d 

dZ l t {u) = Y D i{Zl{")W{t) 

i=0 

reaches x at time t > 0, i.e. Z}(lo) = x. 

Taking for each of the Xj appropriate curves luj with the previous property, we 
obtain - by Taylor expansion of the solution of the non-autonomuous equation as 
in Remark|3|- the familiar version of cubature formulas on Wiener space, namely 

Y e n ---e lk E(f odBi\o...odBil) = 

j eA J0<t!< -<t k <t 

dcg(/)<m 

Y X J E e h'-- e i k dujf(h)---dujj k (t k ). 
j=i ieA J0<ti<-<t k <t 

deg(/)<m 

This leads to a redefinition of cubature formulas, where we replace the points Xi 
by endpoints of evolutions of ordinary differential equations Z\(uji) — Xi for i = 
1 r 

Remark 4. As proposed in [Sj we could also find for any point x € G™ 1 a H 1 - 
trajectory uj : [0,T] — > M. d . such that cu(t) := (t,uj(t)) fort 6 [0, T] satisfies Z\{uj) — 
x, hence we could equally work with this simpler class of curves, which we do not 
do in this article. 
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Theorem 3. Given C 00 -bounded vector fields Vq, . . . , Vd on R , then the diffusion 
process 

dY t y=J2Vi(Yi>)odBi, 

i=0 

yV _ 

for y £ M. N admits the following cubature formula of degree m > 1 : for < t < T 
we find H 1 -curves w\, . . . , tu r : [0, T] — ^ M, d+1 and weights X\, . . . , \ r > 0, such that 

r 

E(f(Y t v )) = £ VWfo)) + 0(t^), 
j=i 

for f £ C?°(M. N ). The constant for the order estimate depends in general on 
derivatives of f up to order m + 1. The curve (Y t v (uj))o<t<T is understood as 
solution of 

d 

dYy{u)=Y J V i {Y?{u))dw i {t), 

z=0 

Remark 5. Notice that in the hypo-elliptic case we can improve Theorem\^ by 
methods proved in 0. Fix m > 1. A family of random variables {Zj for I £ 
A, deg(7) < m} is called m-moment generating if Z$ = 1 and 

E(Z h Z l2 ■■■Z h ) = E{B° h {!)■■■ B o1 ' (1)) 

for all Ij £ A with deg(I\) + • • • + deg(/;) < m. Here we apply 

B o(h {t) = 1 

B oI * {i) {t) := f B oI {s)odB\ 
Jo 

for < t < T. The results on order estimates in 6 apply for m-moment- generating 
families. 

Given now a cubature formula, the coefficients J Q<t < ... <t <1 dio 1 ^ (t\) • • • ckv^ 1 " (tk) 
in the basis • • • ei k , which appear with probability Xj constitute a m-moment gen- 
erating family (see 8 ): this is basically due to the fact that products of iterated 
integrals like B° h (1) • • • B o1 ' (1) can be decomposed as linear combination of iter- 
ated integrals (for instance B\B\ = .B°l 1,2 )(l) + B°( 2,1 \l)), hence the assertion 
on expectations of products follows from the hypothesis on the expectation of it- 
erated integrals. Equivalently one can see this phenomenon also by transforming 
the defining differential equation for the process (Xf) t >o to the Lie algebra in the 
exponential chart. It is remarkable that - due to this observation - we obtain 

r 

E{{Xlf) = Y,HZli"i))\ 
»=i 

for all < k < m, only if the equation holds for k = 1, hence 

r 

E(\og(X 1 1 )) = Y / ^ogZ 1 1 (cu l ) 

i-l 

by definition of the logarithm. 
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3. Calculation of the Greeks 

For the calculation of Greeks we proceed by methods from Malliavin Calculus 
(see ^U] and |§] for all details). We shall consider stochastic differential equations 
of the type 



Y x = y, 

on a stochastic basis (f2, T , P), where we are given a d-dimensional Brownian motion 

(-B|)o<t<r,j=i d in its natural filtration, up to a finite time horizon T > 0. For 

the vector fields V% we shall assume the following regularity assertion: 

• The vector fields 

Vi : R N -> R N 

are C°°-bounded. 
If the distribution 

(3.1) V LA {y) := (Vi(y), . . . , V d (y), [V h Vj}(y), . . . for i, j, ■ ■ ■ = 0, . . . ,d) 

has constant rank N for y E M. N , we say that Hormander's condition holds. Here 
and in the sequel we apply the notion of Lie brackets of two vector fields V, W : 
r n _^ r n^ ^ w ^ ._ dy ^ . W ( y j _ dW(y) ■ V(y) for y E R N . 

Theorem 4. Fix t > 0, v E R N and y E R N . If condition \3.1$ holds, there is a 
weight ir E T)°° , such that for all bounded, measurable functions f : M. N — > R the 
equation 

±\ e=0 E(f(Y t y + n))=E(f(Y t y )*) 
holds true depending on t,v,y and the whole stochastic process (Y' t !/ )o<t<T- 
Proof. For the proof see [H] and JHj- D 

We shall apply usual notions of Malliavin calculus (see [H] and JHI)- The first 
variation Jo^t(y) denotes the derivative of (Y^)o<t<T with respect to y, hence 

d 

dJ ^ t (y) = ]T dVi(Y t v )J Q ^ t (y) ° dB\ 
Jo^o(y) = idN, 

where we apply the Jacobian dVi of the vector field Vi as defined in !2.4l This is an 
almost surely invertible process and we obtain the representation 

D k s Y t y = Jo^ t ( 2 /)(J ^ s (y))-V fe (y/)l [0 , 4] ( s ) 

for < s < T, < t < T, 1 < k < d of the Malliavin derivative. We have the 
fundamental partial integration formula 

d r T 



fe=l 



E( / D k s Y t V a k ds) = E{Y?8{a)), 



n 
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where (a s )o< s <T '■= ( a o< s <T)k=i, ...,d G dom(<5) is a Skorohod integrable process. 
5(a) is real valued random variable. Notice that if a is predictable and square- 
inegrable, then 

d „t 



5(a) = W a k s dB k s , 
i_i Jo 



k=l ' 

hence the Skorohod integral coincides with the Ito integral. 

For the pull-back of the stochastic flow, we can provide a similar stochastic 
Taylor expansion on the space of vector fields as in Section |3 or - in the case of a 
Lie group - on the Lie algebra. We apply the following notation: 



:= V 

■■= [V,Vi\, 

for I = (ii, ...,i k )eA. 

Theorem 5. Given C°° -bounded vector fields Vq, . . . , Vd on K. , then the diffusion 
process 



dYV = J2Vi(Yi>)odBi 

i=0 

yV _ 



fory S M. , admits the following series expansion for the pull back Jo^ t (y) l V(Y^) 
on smooth vector fields V : M. N — > R w . 



J^t(y)- l V(Y t y) = C (D V ^ I odB tl °---°dB. 

Tr=A J 0<tx<-<t h <t 



ll+Rn&V), 



dcg(/)<n 



for < t < T . For the remainder term we obtain 

Rn(t,V) 

E / Jo^ to (y)-^ ,...^m>^ t > 

, n.- „ J0<tn< - <ti.<t 



dog(ii,...,i fc )<n 

(»o,— ,ih)BA 
deg(i ,...,ifc)>n 



for < t < T. 

This Taylor expansion leads on the Lie group G™ 1 to an action on the Lie algebra 
g% v namely 



Xtw(Xl)- x = J2 fei>~ '[ei h M'~] odB\\o-..odBil 
IeA Jo<t 1 <~-<t k <t 

dcg(7)<m— 1 

by inserting the vector fields Di for Vj and evaluating at x — 1 . Notice the decrease 
in order from m to m - 1, since in q™ 1 Lie brackets only up to order m — 1 can 
be non-zero. This stochastic process on the Lie algebra g™ 1 will be applied in the 
sequel in order to construct approximate, universal weights for the calculation of 
the Greeks. First we need a Lemma on scaling invariance. 
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Lemma 1. The following identities in law for the stochastic process (X^)o< s <T 
hold, 

(^^ X s)o<s<l =law (Xg t ) <s<l, 

(a^Kxzmxz)- 1 ]^! = law ((Vt) dc ^»)(x: t )e l (x| t )- 1 ) < s < 1 

for t > and w <G W n (so deg(u>) — n). 

Proof. By scaling invariance we know that (VtBl)o< s <i =i aw (B' l st )o< s <i for t > 
and i = l,...,d. This extends to all iterated stochastic integrals. Notice the 
importance of the degree function, which associates degree 2 to eo, and hence the 
correct scaling property. Notice furthermore that for w € g d n 1 n W n the element 
(Xf)w(X^)~ 1 is an element of the Lie algebra but not necessarily of W n \ □ 

Remark 6. Notice that the stochastic process (-^f)o<t<T is not hypo-elliptic, since 
the direction eg always points into directions where the density does not admit a 
logarithmic derivative. Therefore we only calculate with directions w e 0™i/ ( e o); 
which means directions with vanishing eo component. In all those directions we are 
able to conclude a result on differentiability. 

Remark 7. In the sequel we shall construct Skorohod integrable processes (a*)o< s <T 
by defining them only on [0,t] for < t < T. It is understood that these processes 
are defined on [0, T], but vanish for s > t, hence also a* = a*l[ 0j T]( s ) f or < s < T 
and i = 1, . . . , d. 

Proposition 1. Fix w G Q™ 1 / (eo) (no eo component). For t > there are 
Skorohod-integrable processes (a*)o< s <T (the t-dependence is not visible in our no- 
tation) with a* = a*l[ 0j T](s) f or < s < T, such that 
d 



and 



I (Xj)e i (Xl)- 1 aids = A v - t (w) 
i=i Jo 

j^E{j\<j s fds) = 0{l). 
i=i Jo 



Furthermore we obtain that for any bounded measurable function f : G™ 1 —> K, the 
equation 

±\ e=0 E(f(X: +eA ^ w) )) = E(f(Xn^i) 

holds true, where 7rJ\ := S(s a s ) denotes the Skorohod integral of the strategies 
a. 

Proof. The stochastic process (Xf )o<t<r is not hypo-elliptic, but we prove that 
derivatives in all directions w £ q™ 1 with vanishing e component exist. We can 
define the processes a* by the standard method of Malliavin Calculus: we choose 
a scalar product (,) on which respects the grading, i.e. Wk is orthogonal to 
W\ for k ^ I, and define the reduced covariance operator C* as symmetric operator 

(y.Cy) :=£ f \{Xl)ei{Xl)-\y) 2 ds. 
i=i J ° 
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This operator is almost surely invertible on 0j\/ (eo) with p-integrable inverse (see 
Obviously the image of C* lies in (eo), since there is no eo-component in 
(X])~ 1 eiXl ) the kernel of C* can be calculated by the classical method from [iu] 
and vanishes on Q r ^\/ (eo), hence C'lg^/^o) 1S invertible with p-integrable inverse, 
since the Norris-Lemma applies. 

Furthermore we observe the following scaling property: notice that A^ is a 
self-adjoint operator on g™ 1 for t > 0. We can compare A^C 1 A^ with C* in law 
via Lemma n 

(y,C t y) = J2 f \iXDeiiXD-\y) 2 ds 
i=i ■ /o 

i=i J ° 

c£ ~1 2 
i=l ^° 

= E f {{XDeiiXD-^yf du. 
i=i - 70 

So we have that 

((^ s 1 )e. i (X|))o< s < t ,C*) = law ((^A^KJfJJei^Jj-^Jo^t^^A^) 

for i > 0. We define the (non- adapted) strategies via 

o; := ((XDe^XD.XCX^/^A^w)) 
for 0<s<t, i = 1, . . . ,d, and obtain 

£ f ((Xl)e i (Xl)-\v)alds = 
i=i Jo 

f {{Xl)ei{Xl)-\v) ({Xl) ei {Xl)-\ (Cr^H) ds = 

(«,C*(C t )- 1 A v? ( U ;)>=( U ,A V j( U ;)) 

for u € jj™ r Hence the strategies satisfy the desired equation by duality on q™ v 
Concerning the order estimate for the strategies we observe that 

<(X s 1 )e l (X 5 1 )- 1 ,(C t )- 1 A v ^H) 

^aw^A^KX^e^X^-^A^^^^A^A^H^ 
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so consequently 



= E i? (/ 4 (^( x f) e ^ x fr 1 '( cl ) _lu; ) 2<d (7) 

= X>( /' ((Xl)e t {Xl)-\ (C 1 )-^) 2 du 
i=i Jo 



□ 



Remark 8. The strategies are universal in the sense that they only depend on the 
Brownian motion, the choice of a scalar product on and - in a linear way - 
on the "input" direction w. We denote the Skorohod integrals of those strategies by 
n 1 J L 1 :— 8(a) and call them universal weights. Notice also that E^tt™^ 2 ) = 0(1), 
which is seen by a similar calculation as the previous one. Hence there is a universal 
time dependence, which is achieved by calculating the strategies with respect to the 
direction A^w instead of w. We omit the dependence on w in our notation. 

Example 2. We consider the simplest non-trivial example: m = 2 and d = 2. 
Hence the nilpotent algebra is generated by 

W = (1) , 

Wi = (ei,e 2 ) , 

W 2 = (eie 2 ,e 2 ei,e ) , 



and the Lie algebra by ei,e2,eo, [ei,e2]. The dimensions are 6 and 4 respectively. 
We can solve explicitly 

2 
i=0 

by 

Xf = x + xe 1 Bl + xe 2 B\ + xe a t + xe x e 2 \ j odB\ 1 odB% 2 + 
+ xe 2Sl J J 2 odBl o dB] 2 + \xe\(B\f + \xe\(B\f . 

The pull-back yields a shorter expression, namely 

(XDe^Xl)' 1 = e 1 + [e 2 ,e 1 ]B 2 , 
{Xl)e 2 {Xl)- 1 =e 2 + [e u e 2 ]Bl 
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for < t < T , which can also be calculated directly. Consequently the above equation 
reduces to find a Skorohod integrable strategy such that 

aids = y/iwi 



a 2 ds = Vt' 



w 2 



f (alB 2 s ~a 2 s Bl)ds = tw 3 
Jo 



for w = W\e\ + W2C2 + w-$\e\, ei\- This can be done by introducing a scalar product 
on such that ei, [ei, &2\, eo is an orthonormal basis. We then obtain the 
symmetric matrix 



C l = 



I 



\ 



t 






f Q Blds 
t -J* Bids 

$1 Bids ^{Blf + {Blf)ds 














/ 



with respect to the given basis. So C* is almost surely invertible on {e\,C2, [&i,e2\) 
- as claimed in Proposition^ indeed the determinant is given through 



det(C% lt 



2])) 



{(BlY + {BiY)ds-t{ 



I Bids) 2 - t([ Bids) 2 , 
Jo Jo 



which is positive if and only if one of the two Brownian motions is not vanishing 
identically on [0,t]. The order assertion in t is also nicely visible. 

We consider the first derivative of the function y ^ E(f(Y f y )) for t > 0. By 
Hormander's theorem we know that this function is smooth for all bounded mea- 
surable functions / : R N — > R (see for instance |H] or ^Hl)- We generalize this 
assertion on the one hand, since we leave away Hormander's condition. On the 
other hand we only prove an approximative result with a certain time asymptotics. 

Theorem 6. Fix t > 0, v G and y S R . Furthermore we fix an order m > 1 of 
approximation. We assume that the vector v is a linear combination of Lie brackets 
up to order m— 1 (except the direction Vq), i.e. 

dcg(f) 



E 



r 



I£A\{<£,(0)} 

dcg(/)<m-l 



and define 



E 



dog (J) 

t 2 wi[e ii: [e l2 ,[- 



IGA\{<6,(0)} 
dcg(/)<m— 1 



Then there is universal weight ~n - 
all C 00 - bounded functions f : M. N 



tt™± G T>°° associated to A^q(w), such that for 
-> R the equation 



holds true, where the constant in the order estimate only depends on the first de- 
rivative of f. 
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Proof. Assume / is Ch then 



-\ e=0 E(f(Y t y +ev )) = E{df{Y?) ■ J^ t {y) ■ v), 

where df denotes the differential of the function /. Instead of solving the equation 

d ,t 



u-i Jo 



fe=i 

we take the universal weight iz d n 1 = S(s i— > a s ) of Proposition ^ an d solve the 
equation 

v = E fiJ^siv)-^ (37) - Rt(s, V k ))a k s ds 
k=i Ja 

by the above universal construction. More precisely, we choose real numbers wi for 
I e A\ {0} such that 

v= E ^w/^.^fa.h- ,V ik ]---](y) 



7e.A\{0,(O)} 

dcg(/)<m— 1 



and define 



E * * '^[eitJe^J--- ,e ih ]---]. 



ieA\{0,(o)} 

dcg(7)<m— 1 



Then we take the universal weight 7r := associated to iu and solve consequently 
the above equation. This leads finally to 

±\ t=0 E(f(Yy +ev )) = E(df(Y t v ) ■ Jo-*(tf) • w) 

d ,t 



E(df{Y t y ) ■ J ^ t (y) ■ E / Jo^siyy'VkiY^ds) + 0(t^) 
k=i Jo 

E(J2 fdf{Y?) ■ J^ t {y) ■ J^.{ y )- x V k {YJ/)a\d8) + 0{f^) 

= E(J2 f D k s {Y?)a k s d S )+0{t^) 
Jo 



= E(f(Y t y)n)+0(t^), 

hence the result, since due to Proposition we obtain - by the Cauchy-Schwartz 
inequality and integration with respect to t - an order estimate of the type 0(t ~ ): 
indeed, the remainder satisfies 

£(E f df{Y t y ) ■ J ^ t (y) ■ Ri(s, V k )a h s ds) 
k=i Jo 



< (E#( [\df(Y t y)-J ^ t (y)-R h m ( Sl V k )fd S )HJ2E( [\al) 2 d S ))i 

k=i Jo i=i Jo 

< E M (/ sinds ) h = o(t^). 

u-i Jo 



CALCULATING THE GREEKS BY CUBATURE FORMULAS 



17 



□ 



4. CUBATURE FORMULAS FOR GREEKS 



First we have to get expressions for derivatives of heat equation evolutions on 
free, nilpotent Lie groups, which is an easy task for linear functions by the consid- 
erations of Section [3 



Proposition 2. Fix a linear function A on A™ x , then 

±\ £=0 E(\(X? + ™°)) = X(w x exp(i( eo + i Y. 

i=i 

for w x € x{q 7 ^ 1 / (eo)). Consequently 



(4-1) E(X^) = A^M exp(t( eo + \ J2 e?)) 

i=l 

forw £ 9™i/ ( e o)- 

Proof. We know that the function x i— ► E{\{Xf)) solves the associated heat equa- 
tion. The solution of the heat equation on G™ 1 with initial value A is given by 

d 



A(xexp(t(e + -^e I 2 ))), 
i=i 

hence the derivative in direction w x can be calculated and yields 

1 d 

exp(£(e + - ^ e ?)))' 

i=l 

wherefrom the result follows by duality and Theorem [SI □ 

Next we apply Tchakaloff's theorem twice in order to obtain the appropriate 
cubature result. 

Theorem 7. Fix a free, nilpotent Lie group G™ 1 and w € 0™i/ ( e o), then there 
are points xi, . . . , x r £ G 1 ^ 1 and weights /xi, . . . ,/i r ^ such that 



^ d r 

A^iw) exp(i(e + - ^ e?)) = ]T 



i=l i=i 

Furthermore r < 2 diniR . 

Proof. We write 7r = 7r5\. We define two probability measures on Wiener space, 
absolutely continuous with respect to Wiener measure P, 

dQ+ 1 
dP ~ E{ix + f +1 
dQ- 1 

if the respective positive and negative parts have non-vanishing expectation. Then 
E P (X^) = E(ir + )E Q+ (X?) - E(n_)E Q _(Xf). 
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There is a null set N on Wiener space, such that on N c the process X£ G G d n 1 , hence 
also Xf 6 G d n 1 almost surely with respect to Q±. Consequently by Tchakaloff's 
theorem we find points xf , . . . ,x^ and weights such that 

i=i 

which yields the desired formula by Proposition □ 

Remark 9. We can additionally require that 5ZI=i M« = 0> X)i=i iMil ^ 2 and 
\(H\ <lfori = l,...,r. 

Hence we can formulate a Cubature result for the calculation of Greeks. 

Definition 2. Fix m > 1 arjrf < t < T, a cubature formula for the first derivative 
in direction A^/jw £ 0™i/ (eo) is given by a finite number of points x±, . . . ,x r € G^i 
and finitely many weights fi±, . . . , fi r ^ 0, such that 

r 
3 = 1 

or equivalently due to formula j4-l\ ) 

^ d r 

A vtH ex P^( e o + 2 ^Z e D) = ^2^- 

1=1 j=l 

Again we use trajectories Ui to represent the points Xi as endpoints of evolutions 
of ordinary differential equations Z t {uJi) = Xi for i = 1, . . . , r. Furthermore r < 
2dim R A^ 1 . 

Theorem 8. Given C 00 -bounded vector fields Vo, . . ., Vd on Mr , then the diffusion 
process 

d 

dY? = Y. vl{ yt v ) odB i 

for y E M. N admits the following cubature formulas for derivatives in direction 
v S R^. Fix m > 1 and assume that 

(4-2) v= Yl t^wrlVi^ViA--- ,V ik ]---](y), 

/e.A\(0,(o)) 

dcg(7)<m— 1 

then we obtain 

l\ e=0 E(f(Y t y + n) = j^nfiy?^)) + o{t^i 

3=1 

taking a cubature formula, which was derived in G™ 1 and C°° -bounded f . The 
constant depends in general on derivatives of f up to order m+l. We can determine 
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the weights fij and the trajectories Wj by 

1 d 

A^Hexp(i(eo + -^ e 2)) 

i=l 

= y]Mi e n---e lk dw^(tt)---dujf(t k ), 



i=i ieA\($,(o)) 

dcg(/)<m 

/£i\(»,(o)) 

dcg(/)<m— 1 

Remark 10. Since we shall in practise apply this procedure to smooth functions 
- as discussed in the recipe after formulas \ and \1.J$ - this already yields the 
interesting result for applications. 

Proof. By the previous constructions we obtain 

e{x^2i) = X> E • • • e H / ^ (*0 • • • ^? (**)• 

i=i i e A Jo<t 1 <---<t k <t 

deg(/)<m 

We then insert instead of the vector fields Vi and apply those vector fields to the 
function / at y € R . Then we obtain on the left hand side - due to stochastic 
Taylor expansion - and on the right hand side due to Taylor expansion of the 
non-autonomous equation, 

r 
i=l 

The left hand side yields then the order estimate by partial integration and Propo- 
sition the right hand side by the fact that iterated integrals with respect to ujj 
behave like t~^~ . □ 



We additionally have an assertion on the construction of cubature paths due to 
the scaling properties. Assume that we found trajectories LOi and weights ^ for 
m > 1, d > 1 and t = 1 fixed, then we can construct solutions for all t > with 
the same m, d: the equation 



1 d - f 

( e o+oE e ')"E^ E e *i"' e H / duf (h) ■ ■ ■ dajf (t k ) 

1 i=i i=i ieA ^o<t!<-<t fc <i 



w exp 



holds, hence the trajectories 



re 

dcg(/)<m 



•#(st) :=^°( S ), 
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for i = 1, . . . , d, j = 1, . . . , r and < s < 1 satisfy 

1 d 

A^Hexp(i(e + -^e2)) 



2 

i=l 



i=l 16.4 70<ti<-<t, 



<t 



dcg(/)<m 



since 



y^Mj J! e n---e lk drjf (h) ■ ■ ■ dr)f (t k ) 

i=i i e a Jo<ti<-<t k <t 

deg(7)<m 

i=l JC.4 iO<t!<-<t fc <l 



/e.4 

dcg(/)<m 



^ J0<t, < — <th<l 



dcg(/)<m 

1 - 

= A v7 (w exp(e + - X e ? )) 



1 d 

A^H e X p(t(eo + -^e?)) 



2 ■ 
2 



i=l 

This leads to the desired assertion. 

5. Applications 

5.1. Example for m = 1 and d > 1. In [1] the authors provide the following 
expression for the Malliavin weight ir: given C°°-bounded vector fields V\, ■ ■ . , Vd 
on M. d such that a uniform ellipticity condition holds, i.e. there is 6 > such that 

d 
i=l 

for all £ € R d , with respect to a standard scalar product on R d , the formula 
(5-1) || e=0 £tf(y/ + "')) = £(/0if)i j\a-\Y^){J^ s {y) ■ v)) T dB s ), 

where we understand dB s as column vector, where the random matrix a~ 1 (Y^) 
acts on. a is defined via 

<r(v) ■= {Vi(y),...,V d (y)). 

Notice that this formula provides adapted strategies, therefore the Skorohod in- 
tegral can be replaced by the Ito integral. We shall compare this formula to the 
formula obtained in our setting. We choose m = 1, we can calculate ir^ 1 explicitly, 
namely 

d 

* 1_ § ^ 
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for i = 1, . . . , d and w has ei-component w% for i = 1, . . . , d. Hence we approximate 

d 



^-\ e=0 E(f(Yr eVlv )) = E(f(Y t y ) £ bQ) + 0(t), 

i—l * 

where stem from the solution of the equation 

d 
i=l 

consequently w = a^ 1 (y) ■ v. Hence 

±\ e=oE (f(Y t y + ^)) = E(f( Y y) J2 B f~ 1{ j}' v ) + °®> 

i—l v 

lu =0 E(f(Ym) = E(f(Y?) £ B f~ 1{ f' V ) + 0(Vi 

i—l 

which is obviously the expansion of the precise formula l|5.1|) , since 

,t d 



E{f{Y?)\ f (a-HY^(Jo^(y)-v)) T dB s )=E(f(Y t y )^Bi ^^ )+ 

1 JO 1 



= 1 

t 



+E(f(y?)i I (cT- 1 (Yy)(J ^ s (y)-v)-a- 1 (y)-v) T dB s ) 



i=l 

t 



for / a C^°-function. 



+E{\ I D s f(Yy)(a- 1 (Yy)(J ^ s (y)-v)-a' 1 (y)-v) T ds) 



E(f(Y t ^Bt a 1{ f' V ) + 0(y/i) 

i=l 



5.2. Example for m = 2 and d = 2. In order to demonstrate the method we 
calculate one example of the above methodology. First we fix m = d = 2, hence we 
can include the calculations of Example [2 We first fix a direction w = e%. Then 
we find one solution of 

1 2 

Vtei = Vtei exp(t(e + - £ e 2 )) = £ /i,-^, 

i=l i=l 

which can be given through r = 2, fx± = ^, /i 2 = ~\ and X\ = exp(\/tei), X2 = 
exp(— vtei). Consequently the trajectories are given by 



w i( s ) = I ft 




W 2<» = ( 
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which leads to 

Zl(ui) = 1 + Vie! + = exp(Vtei), 

Zl(w 2 ) = 1 - Vie\ + ^-e\ = cxp(-Viei). 

This in turn leads to the following approximation scheme: given vector fields 
Vo,Vi,V2 on M. N (notice that Hormander's hypo-ellipticity condition (|1.2|) is not 
necessarily fulfilled): 

• solve the ordinary differential equations 

which yields flows in the diffusion direction V\ starting at y £ M. N . 

• fix the direction for the directional derivative due to the above construction 
e\ corresponds via formula (|4.2|l to V\, hence v = ViV\{y). 

• the derivative can be approximated through 

±\ e=0 E(f(Y t y +£VIVl{y) )) = \{m v {ux)) - f (??(«»))) + o{th 

where the constant in the order estimate depends on the third derivative 
of/. 
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